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Using equations, governing np pairing correlations in S = 1, T = pairing channel 
(PRC 63 (2001) 021304(R)), it is shown that at low densities equations for the energy gap 
in the spectrum of quasiparticles and chemical potentials of protons and neutrons allow 
solutions with negative chemical potential. This corresponds to appearance of Bose- 
Einstein condensate (BEC) of deuterons in low density region of nuclear matter. 



The transition from BCS superconductivity to Bose-Einstein condensation occurs in a 
Fermi system, if either density is decreased or the attractive interaction between fermions 
is increased sufficiently Recently it was realized that this phase transition takes place in 
symmetric nuclear matter, when np Cooper pairs at higher densities go over to BEC of 
deuterons at lower densities [ EJ 12] • During the phase transition the chemical potential 
changes its sign at certain critical density (Mott transition), approaching half of the 
deuteron binding energy at ultra low densities. Here for studying corresponding phase 
transition in asymmetric nuclear matter we shall use equations, obtained in Ref. for 
description of np pairing correlations in S — 1,T — pairing channel: 

A(k) = -knk,k')^(l-/(^)-/« (1) 

9 = fEO-g ^-m)-f^)])^^n k , (2) 

<*Q = | E(/(^ + ) - /(£*)). E t = E k ± ^03 = y/? k + A*(k) ± /xoa, (3) 

where f(E) is Fermi distribution, e k is the single-particle spectrum, /xoo,/Uo3 being half 
of a sum and half of a difference of proton and neutron chemical potentials. Eq. (|TJ 
is equation for the energy gap A and Eqs. (j2J, © are equations for the total density 
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Q — Qp + Qn and neutron excess 5g = g n — g p = ccp (a being the asymmetry parameter). 
Then, introducing the anomalous density 

m =< >= - ) - /to) 

and using Eq. (J2J), one can represent Eq. (pQ) for the energy gap in the form 

-^(k) + (1 - n k ) £ V(K k')V(k') = 2/i 00 ^(k). (4) 
m V 

In the limit of vanishing density, n k — > 0, Eq. (j^J) goes over into the Schrodinger equation 
for the deuteron bound state [13 El- Corresponding energy eigenvalue is equal to 2/j, 00 . 

Further for numerical calculations we shall use the effective zero range force, developed 
in Ref. [[Sj to reproduce the pairing gap in S = 1, T = pairing channel with Paris NN 
potential: 

V{n, r 2 ) = ^o{l - rjf^i^y^ - r 2 ), (5) 

where go is the nuclear saturation density (go = 0.16 fm~ 3 ), Vo,r),P are some adjustable 
parameters. Besides, in the gap equation (JTJ), Eq. (JHJ) must be supplemented with a cut-off 
parameter e c . We utilize the following set of parameters: r\ = 0, v q = —530 MeV-fm 3 , m = 
m G, e c = 60 MeV, where mo is the effective mass, corresponding to the Gogny force D1S [ 




Figure 1. Chemical potential /zoo as a function of: (a) density at zero temperature, 
(b) temperature at fixed densities (low density region). 



E]. This set was developed to produce the bound state between two extreme values: at 
zero energy and at the deuteron binding energy [UJ. 
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Figure 2. Energy gap as a function of density at zero temperature and different asymme- 
tries. 



First we consider the case of symmetric nuclear matter (a = 0). In Fig. 1(a) it is shown 
the zero temperature behavior of the chemical potential /ioo- It is seen that the chemical 
potential at some density gt (@b « 3-1CT 4 fm~ 3 ) changes its sign, that, according to Eq. (JH), 
corresponds to appearance of deuteron-like bound states in nuclear matter. When density 
tends to zero, the chemical potential approaches its asymptotic value /ioo = ~£fe/2, being 
the binding energy of a bound state. Therefore, we can conclude that under lowering 
density np superfluidity smoothly evolves into the BEC of np bound states (deuterons). 

In Fig. 1(b) we present the results of numerical determination of the temperature de- 
pendence of the chemical potential /i o for the fixed values of density. If density is low 
enough then under decrease of temperature the chemical potential becomes negative (the 
curve 2), that corresponds to formation of bound states. At very low densities np con- 
densate exists only in the form of BEC of deuteron-like bound states (the curve 3). If 
densities is, however, high enough, np Cooper pairs survive even at zero temperature (the 
curve 1). 

Now we consider asymmetric nuclear matter. The results of numerical calculations for 
the energy gap as a function of density for different asymmetries at zero temperature are 
shown in Fig. |21 As one can see, with increasing asymmetry the value of the energy gap 
decreases and the density interval, where np pairs exist, shrinks in the direction of zero 
density. In reality solutions exist for any a < 1 (the phase curves for larger values of a 
are not shown in Fig. |2J) and corresponding density interval contracts more and more to 
the point g = 0, when asymmetry increases. To understand why the isospin asymmetry 
loses its efficiency in destroying np pairing correlations in low density region, let us note 
that at zero temperature the contribution to the integral in the gap equation gives the 
interval [0, e c ], excluding the domain [/i o — Ae, /too + Ae], where Ae = J Uq 3 — A 2 . In the 
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Figure 3. Various chemical potentials as functions of density at zero temperature and 
different asymmetries: (a) //qo, A*03, Hp, Hm (b) Hp (lower curves), //„ (upper curves). 



weak coupling regime (/xqo ^ A) with increasing asymmetry the width of this domain 
also increases, that results in considerable reduction of the energy gap magnitude, until 
it completely vanishes. However, when the chemical potential passes through zero and 
becomes negative, only part of the window participates in suppressing the energy gap, 
with the right end of the blocking interval going to zero at g — > 0. 

In Fig.|2fa) it is shown zero temperature behavior of chemical potentials /ioo, Has, Hp, Hn 
as functions of density at very low densities of nuclear matter and finite isospin asymmetry. 
One can see, that at some critical density the chemical potential noo changes its sign and 
np Cooper pairs smoothly go over into deuteron bound states. The asymptotic behavior 
of chemical potentials at g — > is Hoo, Hos ~^ — £ b/2, and, hence, /i p — > — e&, n n — ► 0. 
This asymptotic behavior does not depend on isospin asymmetry, that is confirmed by 
the results of numerical calculations, shown in Fig. 0Kb), where the density dependence of 
chemical potentials Hp, Hn is depicted for different values of the asymmetry parameter. 
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